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Abstract 

X \ The radiative decay K~ -^ /i~i^y^7 has been studied at ISTRA+ setup in a new kinematic 

I ■ region. About 22K events of K~ -^ /^"i^/^T have been observed. The sign and value of 

^ Fv — Fa have been measured for the first time. The result is Fv — Fa = 0.21±0.04(stat)± 

^: QM{syst). 
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1. Introduction 



Radiative kaon decays are dominated by long distance (low energy) physics. For low 
energy processes there are no direct predictions from SM and effective theories such as 
O ' Chiral perturbation theory (xPT) are used. xPT gives quantitative predictions for most 

kaon decay modes. That is why radiative kaon decays provide a testing ground for xPT. 
Moreover, these decays are sensitive to New Physics. 
^ ■ The decay K~ -^ M~^/i7 is sensitive to hadronic weak currents in low-energy region. 

The decay amplitude includes two terms: internal bremsstrahlung (IB) and structure 
dependent term (SD). IB contains radiative corrections from K~ -^ l^~i^i^- SD allows to 
probe electroweak structure of kaon. 

The differential decay rate can be written in terms of standard kinematic variables 
X = 2E^/Mk and y = 2£'^/M/e (see [1] for details), E^ being photon energy and E"^ muon 
energy in cms. It includes IB, SD^ parts and their interference INT^. The SD^ and 
INT^ contributions are determined by two formfactors Fy and Fa- 

The general formula for decay rate is as follows: 

^^ - AiBfiB{x,y) + Asd[{Fv + FAffsD+{x,y) + {Fy - FAffsD-{x,y)] 



dxdy 
Aint[{Fv + FA)fiNTA^,y) + {Fv - FA)fiNT-{x,y)] 

where 
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fsD+{x,y) = [x + y- l-r][{x + y- l)(l-x) -r], 
fsD-ix,y) = [1 - y + r][{l - x){l - y) + r], 
fiNT+{x,y) = [ ^(^ + y-l- r) ][(l ~ ^)(1 -x-y) + r], 

fiNT- (x, y) = [^n^^^^^f:rf—A]i^^ - {l-x){l-x-y) -r], 



and r = 
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AiNT = ^Kf,2'T~T\ — n2 ~f^> I^ these formulae, a is the fine structure constant, F^ is 

/C+ decay constant {Fk = 156.1 ± 0.2 ± 0.8 ± 0.2MeV [2]) and Tk^^ is ^^2 decay width. 
Dahtz-plot distributions for different terms are shown in Figs. [U-Ill 

Fv ± Fa are calculated within xPT (O(p^) [1], 0{p^) [3]) and LFQM model [4]. In 
general, Fy and Fa depend on g^ = (p^ _ p^)2 = mI{1 - x). In the 0(p'^) xPT they 
are constant and Fy + Fa = 0.137, Fy — Fa = 0.052. We will initially assume Fy and 
Fa constant and then test for their dependence on q"^. 
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Figure 1: IB. 



Figure 2: INT". 
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Figure 3: SD+. 



Figure 4: SD". 



Experimentally, the decay K -^ l^^iil was studied mostly in the IB dominated 
region (see [5(],J6|,J7|). There was only one formfactor measurement in E787 exper- 



iment |8|]. In this study, SD+ term was extracted and \Fy + Fa\ was obtained to 
be \Fv + Fa\ = 0.165 ± 0m7{stat) ± Omi{syst). Also Fy - Fa was constrained: 
—0.04 < Fy — Fa < 0.24. Fy — Fa was measured by E865 experiment in K ^ fiiye^e~ 
decay [9]: Fy — Fa = 0.077 ± 0.028. The goal of our study is to measure K~ -^ l^~^iil 
decay in the kinematic region where INT~ term (and hence Fy — Fa) can be extracted. 



2. ISTRA+ experiment 

2.1. Experimental setup 
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Figure 5: Elevation view of the ISTRA+ detector. 



The experiment was performed at the IHEP 70 GeV proton synchrotron U-70. The exper- 
imental setup ISTRA+ (Fig.[5j) was described in details in llO|]. The setup was located in 
the negative unseparated secondary beam. The beam momentum in the measurements 
was ^ 26 GeV/c with Ap/p ^ 1.5%. The admixture of K~ in the beam was ^ 3%. 
The beam intensity was ^ 3 • 10^ per 1.9 sec U-70 spill. The track of a beam particle 
deflected by Mi was measured by BPCi — BPC4 (1mm step multiwire chambers), the 
kaon identiflcation was done by Cq — C2 threshold (7-counters. A 9 meter long vacuum 
decay volume was surrounded by 8 lead glass rings LGi — LGs used to veto low energy 
photons. SP2 was a lead glass calorimeter to detect/veto large angle photons. Tracks of 
decay products deflected in M2 with ITm fleld integral were measured by PCi — PCs 
(2mm step proportional chambers); DCi — DCs (1cm cell drift chambers) and flnally by 
2cm diameter drift tubes DTi — DT4. Wide aperture threshold Cerenkov counters Ca, 
C4 were flUed with He and were not used in the measurements. Nevertheless Cs was used 
as an extension of the decay volume. SPi (ECAL) was a 576-cell lead glass calorimeter, 
followed by HC (HCAL) - a scintillator-iron sampling hadron calorimeter. HC was sub- 
divided into 7 longitudinal sections 7x7 cells each. MH was a 11x11 cell scintillating 
hodoscope used to improve the time resolution of the tracking system, MuH was a 7x7 
cell muon hodoscope. 

The trigger was provided by Si — Ss^ S^ scintillation counters, Cq — C2 Cerenkov 
counters, analog sum of amplitudes from the last dinodes of the SPi : Tq = Si • S2 • Ss • 
Co • Ci • C2 • S^ • 2(5^ Pi), here S^ was a counter downstream the setup at the beam focus; 
S(S'Pi)- a requirement for the analog sum of ECAL amplitudes to be above ^3 GeV. 
The last requirement served to suppress the K -^ /xz/_ decay. About 10% events were 
recorded with a different trigger: Ti = Si • S2 • Ss • Co • Ci • C2 • S^. This prescaled trigger 
allowed to calculate trigger efl&ciency as a function of the energy released in ECAL. 



2.2. Data and MC samples 

We use high-statistics data coUected in Winter 2001 run. About 332M events were 
stored on tapes. This statistics was complemented by 200M MC events generated with 
Geant3 [ul]. The MC generation included realistic description of all ISTRA+ detectors. 

3. Event selection 

3.1. Selection criteria and general cuts 

The decay signature is defined as follows: one primary track (kaon), one negatively 
charged secondary track identified as muon; one shower in ECAL not associated with the 
charged track. Muon identification using ECAL and HCAL is described in our previous 
papers (H^H). 

Several cuts are applied to clean the data: 

• number of beam and decay tracks in both X and Y projections is equal to 1; 

• CL (track quality) of primary tracks in both X and Y projections must be greater 
than 10"^; 

• CL of decay tracks is greater than 0.1 (decay-X) and 0.15 (decay-F); 

• the angle between primary (kaon) and secondary (muon) track is greater than 2 
mrad. 

The last cut eliminates most of undecayed beam particles. The quality of decay track 
(described quantitatively by CL) is worse than that of beam track because of multiple 
scattering and detector resolution. 

Cuts containing photon energy include: 

• gamma energy in kaon rest frame is greater than 10 MeV; 

• no photons in SP2 calorimeter (energy threshold is 0.5 GeV for total energy release); 

• no photons in GS. 

For vertex characteristics we have the following requirements: 

• z-coordinate must be within the interval 400 < z^tx < 1600cm; 

• (-3) < Xytx < 3cm; 

• (-2) < y^tx < 6cm; 

• CL of general vertex fit is greater than 10~^. 
Additional cuts are applied to suppress backgrounds: 

• number of hits in matrix hodoscope (MH) is less than 3; 

• missing momentum p ^^55 =Pk— P/i— Pj does not point to the ECAL central 
hole (this cut eflFectively rejects background from K -^ nn^ decay when the lost 
photon from n^ -^ 77 goes into the hole). 



3.2. Trigger efficiency 

As To trigger described in Section 2 contains energy threshold in SPi the trigger 
efficiency as a function of energy released in ECAL could be found using events with Ti 
trigger: Ctrg = (Ti H^o) / Ti. Trigger curve is shown in the Fig. [61 The fit is done using 
Fermi function. For the further analysis only events with Tq are kept and these events 
are weighted by the factor of l/ctrg- 
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Figure 6: To trigger efficiency. Points - data, curve - fit by Fermi function. 



4. Signal extraction 

Distribution over M(/xz/7) is used for signal observation. M^(/xz/7) = (P^ -\- Py-\- P^Y 
where P^^P^^ P^ are 4-momenta of corresponding particles; missing mass m^ is supposed 
to be equal to so that p ly = p k — P ^ — P-i]Ey = \piy\. M{iiiiyj) peaks at K~ 
mass for the signal. Main background comes from 2 decay modes: K~ -^ ii~ vn^ {K jiZ) 
and K~ -^ ti~ti^{Kti2) with one gamma lost from tt^ -^ 77 and tt misidentified as ji. 
Dalitz-plot distributions for X/iS and Kti2 are shown in Figs. [71 O 





Figure 7: Dalitz-plot density for i^/i3 bkg. 



Figure 8: Dalitz-plot density for i^7r2 bkg. 



^.7. Signal extraction procedure 

The procedure starts with dividing aU kinematic (x^y) region into stripes on x (x- 
stripes). The stripe width is Ax=0.05 (AE"* ^ 24:MeV). In every x-stripe we put a 
cut on y: y^^in < y < ymax- Vmin ^ud ymax ^rc selected from the maximization of signal 
significance defined as -mr^- 

Besides distributions over M{iJiv^) and y, we use cos 0*^ for the signal observation, 
0*^ being the angle between p^ and p^ in c.m.s. We put a cut on cos 9^^ to reject 
background in those stripes where distributions over cos 0*^ for signal and background 
differ a lot (for example, in Fig. [13] in the the distribution over cos 9^^ signal peaks at 
0.8 and background peaks at -0.4). 

Finally for each x-stripe we obtain events with cuts on y and cos 0*^. Now we 
construct M(/xz/7) which will be used for the fit. Fitting M(/iz/7) alone is not sufficient 
because in some stripes distributions for signal and background are very similar. Also it 
would be difficult to distinguish between two backgrounds - KfiS and K7r2. That is why 
we take three histograms (y; cos 0*^ with cut on y] M{iLiiy^) with cuts on y and cos 0*^) 
and fit them simultaneously. Both signal and background shapes are taken from MC. 
MC histograms are smoothed and the result is stored as f{z) function {z = M(//z/7), 
y or cos 9'^^). For better fit we allow these functions to be slightly widen and shifted. 
We do it by using f{k^z + b) instead of f{z) in the fit, where fit parameters k and b 
are the same for signal and background and are different for M(/xz/7), y and cos 9~^^. In 
Figs. [91 and [TOl extra width k and extra shift b for y-histogram are shown. For all selected 
x-stripes k~l and b^O, i.e. our MC describes data properly. 





Figure 9: Extra width ky vs x. 



Figure 10: Extra shift by vs x. 



The simultaneous fit gives signal event number in each x-stripe. As we use the same 
data several times we should take care about correct estimation of statistical error: 

• do simultaneous fit of three histograms and obtain {pi} - best parameter values 
(they correspond to global x^ minimum); 

• take {pi} as initial values and perform y^ /n.d.f. and error estimation for one his- 
togram M(/xz/7) using MINOS program (l4l |. 

4^2. Selected kinematic region 

For further analysis we have selected eleven x-stripes in the following region: 0.05 < 
X < 0.6 {12MeV < E^ < lASMeV). The twelfth x-stripe is used for systematics study 
only. Cuts on y and cos 9^^ are summarized in tablelH 



6 



strip 


cut on X 


cut on y 


Ay 


cut 


on cos 6**^ 


01 


0.05 < a; < 0.1 


0.9-1.1 


0.2 


> - 


-0.8 


02 


0.1 < a; < 0.15 


0.9-1.1 


0.2 


> - 


-0.8 


03 


0.15 < a; < 0.2 


0.85-1. 


0.15 


> - 


-0.8 


04 


0.2 < a; < 0.25 


0.8-0.95 


0.15 


> - 


-0.2 


05 


0.25 < a; < 0.3 


0.75-0.9 


0.15 


> - 


-0.3 


06 


0.3 < a; < 0.35 


0.72-0.87 


0.15 


> - 


-0.4 


07 


0.35 < a; < 0.4 


0.65-0.85 


0.2 


> - 


-0.3 


08 


0.4 < a; < 0.45 


0.62-0.85 


0.23 


> - 


-0.5 


09 


0.45 < a; < 0.5 


0.57-0.8 


0.23 


> - 


-0.7 


10 


0.5 < a; < 0.55 


0.52-0.75 


0.23 


— 




11 


0.55 < a; < 0.6 


0.48-0.7 


0.22 


— 




12 


0.6 < a; < 0.65 


0.45-0.65 


0.2 


- 





Table 1: Cuts on y and cos 0* in x-stripes. 



The ?/-width changes from stripe to stripe, in average A^ ^ 0.2. Our kinematic region 
is sensitive to INT~ term (Fig. [TTl ) and complementary to that of previous experiments 
[8|],|5| (Fig. [12]). Stripe borders are slightly out of allowed kinematic region because of 
resolution. 



X 




Figure 11: INT dalitz-plot density and se- 
lected kinematic region. 




Figure 12: ISTRA+ (green); BNL E787(red 
hatch); KEK-104(blue hatch). 



Results of simultaneous fit for stripes #4 (0.2 < x < 0.25) and #9 (0.45 < x < 0.5) 
are shown in Figs. [131 El The total number of unweighted signal events is ^22K. 



5. Fy — Fa measurement 

For each x-stripe we have experimental event number A^exp from fitting the data and 
IB event number Njb from MC (see Fig. [15]). Then we plot N^xp/NjE as a function of x 
where each bin corresponds to a certain x-stripe (see Fig. [T6l) . 

For IB only we would have Nexp/Njs ^1. It is the case for small x where IB is 
dominated and INT~ is negligible. For large x we see that N^xp also contains negative 
interference term. We fit N^xp/NjE distribution with {fiB{x) — fiNT-i^-iP))/ fiB{x) where 
fiB and fiNT- give MC event number for a certain x-stripe and fit parameter p equals 
to Fv — Fa {Fy and Fa are initially assumed to be constant). The result of the final fit 
is as follows: Fy — Fa = 0.21 ± 0.04(5tat). Number of 'missing events' due to negative 
INT~ term is 1483 which is ^ 3% of expected IB contribution (49722 weighted events). 
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Figure 13: Simultaneous fit in stripe 4: M^jiv^)^ y and cos ^* . Points with errors - data, blue - i^/i3, 
red - i^7r2, green - signal, red line - fit (signal+background) . x^ /^-d-f -=^^7.3/91 (for mass histogram 
only, see text). 
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Figure 14: Simultaneous fit in stripe 9: M(/iz/7), y and cos 0^^. Points with errors - data, blue - i^/i3, 
red - K7r2, green - signal, red line - fit (signal+background). x^/^-<^-/-=141.7/91 (for mass histogram 
only, see text). 



6. Systematic error estimation 

The main potential sources of systematic error are: 

• possible non-ideal description of signal/background shape in the simultaneous fit; 

• cut on X (i.e. number of selected x-stripes); 

• x-binning (i.e. x-stripe width); 

• cut on y in x-stripes; 

• cut on z-coordinate of the vertex; 

• possible contribution of INT+ term. 

Each source is investigated separately and errors are considered to be independent. 
Possible non-ideal description of signal/background shape in the simulta- 
neous fit. For estimation of shape systematics we scale errors in each bin of the final fit 
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Figure 15: x-spectrum. Red points - data, blue 
histogram - IB. 



Nexp/NiB^"^ 




Figure 16: ^exp/NiB for x-stripes and final fit. 

xV^.^./.=7.8/9. 



of Fig. [T6l proportionally to yx^ (x^ is obtained from simultaneous fit in a bin). Then we 
repeat final fit. New value of Fy — Fa is consistent with the main one and the fit error is 

larger: cffu -- 5 x 10"^ We treat cffu as follows: cffu = J (jltat + ^Lst fit ^i^h (Jsystjit be- 



systjit 

ing systematical error caused by non-ideal shape of signal and background distributions: 

^systjit ^ 3 X 10"^. 

Cut on X. Each x-stripe has the width Ax=0.05. By adding/removing stripes 
involved in the fit on the left (right) border and repeating final fit we can see how Fy — Fa 
depends on the x-cut value. For the left border, we take results of 3 fits which include 
stripes 1 — 11 (main fit), 2—11 and 3—11. For the right border, we choose fits including 
stripes 1 — 10, 1 — ll(main fit), 1 — 12. 
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Figure 17: Systematics of cut on x. Left bor- 
der. 
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Figure 18: Systematics of cut on x. Right bor- 
der. 



The resulting plots Fy — Fa vs x-cut are shown in Figs. [TTl and [T8l For the conservative 
estimate of systematics we fit these plots with straight lines. The line slope multiplied 
by the resolution in x (which is taken from MC) gives systematic error of this cut. The 
systematic error of the right border is found to be ^ 1.2 x 10~^ and that of the left border 
is negligible. 

X-binning. Changing width of x-stripes could shift the result of the final fit. To 
test that, we repeat the whole procedure (choosing cuts on y in x-stripes, simultaneous 
fits and a final fit) for two diflFerent values of Ax: Ax=0.035 and Ax=0.07. The smaller 
value of Ax equals resolution in x. By comparing results obtained in the fits with new 
Ax with the main one (Ax=0.05) we find tsyst ^ 2 x 10~^. 
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Cut on y in x-stripes. To investigate this source of systematics we choose cut on 
y in a different way: instead of using significance we select events inside FWHM in y- 
distribution for signal MC. Such cuts on y are stronger than those made using significance. 
We redo simultaneous fit in x-stripes and final fit. The obtained result is consistent with 
the main one. No systematics is found here. 

Cut on z-coordinate of the vertex. To study this systematics we divide events 
into two groups - with z^tx < HOO cm and z^tx > HOO cm. The events with z^tx < HOO 
cm use PCI in the decay track reconstruction while events with z^tx > HOO cm do not. 
Besides that the second group of events has the vertex inside the decay volume filled 
with He. It could be a possible source of systematics. Repeating the whole procedure 
(simultaneous fit in x-stripes and final fit) we obtain two values for Fy — Fa which are 
averaged. The obtained values are compatible. No systematics is found here. 

Possible contribution of INT+ term. Contribution of INT+ term is proportional 
to Fv + Fa' In E949 experiment [8|] only the absolute value of Fy + Fa was measured. 
Fig. [T9l shows that INT+ could be rather large compared to INT~, especially for small 
X. Fortunately in this case INT~ is negligible with respect to IB. Since we know only 
\Fy + Fa\ we add INT+ term with both signs to the fitting function for the final fit and 
treat the difference in obtained value of Fy — F^ as a systematic error. This way we get 

Csyst ^ 1.4 X 10-^ 



N so 




Figure 19: INT + (magenta) and INT~(cian) contributions with Fy + Fa=0.165 and Fy - Fa=0.21. 



Total systematic error. Now we quadratically sum all sources supposing the errors 
to be independent and obtain Cgyst ^ 4 x 10~^. Contributions are summarized in table [21 



source 


systematic error 


non-ideal shape 


3x10-^ 


cut on X 


1.2x10-2 


x-binning 


2x10-2 


cut on y 


- 


cut on Zytx 


- 


INT+ contribution 


1.4x10-2 


total 


4x10-^ 



Table 2: Systematic errors. 
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7. Final result 



With this estimation of systematic error we finaUy get our result: Fy — F^ = 0.21 ± 
0.04:{stat)±0.04:{syst). It is ^ 3a larger than theoretical prediction within xPT at O(p^). 

The 0{p^) xPT gives linear dependence of Fy and Fa on q^ (see [31]) and hence on x. 
We use Fy and Fa parametrization given in [15]: Fy = Fy{0) [1 + A(l — x)], F^=const. 
This theoretical prediction was tested in three ways. First, we take both Fy and Fa 
from 0(/) xPT (Fy(0) = 0.082, F4=0.034, A=0.4) and do the final fit. x^ of this fit is 
21.1/10 (- 2.5a from x^ = !)• Second, we take Fy{0) and Fa from 0(/) xPT and take 
A as a fit parameter. It gives A = 4.0 ±1.0 with x^=8.8/9 (Fig. [20]). And finally we fix 
Fy{0) from 0{p^) xPT and take A and Fa as fit parameters. Correlation between them 
is shown in Fig. [221 Theoretical prediction is slightly out of 3a-ellipse. 

In LFQM, Fy and Fa depend on q^ in a complicated way (see [4]). Final fit is shown 
in Fig. [in LFQM is disfavored (^ 3a from x^ = 1) although can not be excluded. 



Nea:p/NjB^ 




NexplNiB 




Figure 20: xPT 0(/) fit, Fy(0) and Fa taken 
from theory. x^/^-<^-/-=8-8/9. 



Figure 21: LFQM fit, Fy and Fa taken from 
theory. xV^-^-/-=24.1/10. 
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Figure 22: xPT 0{p^) fit, Fy{0) taken from theory x^/n.d.f. = 7.5/8. 



8. Conclusions 

The radiative decay K~ -^ l^~^^l has been studied using in-flight decays at ISTRA+ 
setup. About 22K events of K~ -^ l^~^^l (it is the largest statistics for this decay) have 
been found in a new kinematic region. The negative INT~ term has been observed and 
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as a result Fy — Fa has been measured: Fy — Fa = 0.21 ± 0.04(stat) ± 0.04{syst). The 
result is ~ 3a above 0{p'^) xPT prediction. 

An alternative analysis has been done by our collaboration (preliminary results are 
presented in [16]). The results are in reasonable agreement (x^=1.12). 

Authors would like to thank C.Q. Geng and E. Goudzovsky for the code plotting 
formfactors in LFQM. The work is supported by Russian Foundation for Basic Research 
(grants 10-02-00330 and 08-02-91016). 
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